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The distribution of the fractal dimension of the two-dimensional Ising model at the critical 
temperature measured by the Monte-Carlo simulation is discussed. At small spatio-temporal scales 
it exhibits a multifractal behavior and is well fitted by a non-Gaussian distribution derived from the 
Renyi or Tsallis entropy. 
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Recently a large-deviation analysis Q for the Monte-Carlo data of the two-dimensional Ising model at the critical 
temperature has been attempted. Inspired by this attempt we perform a multifractal analysis for the same model in 
this Short Note. Here we do not distinguish the large-deviation and the multifractal schemes. || 

Multifractal behaviors of the Ising model are expected for random or chaotic lattices, j| while the model on a 
regular lattice is understood in terms of the mono-fractal analysis. H Even for the model on a regular lattice, we 
can expect multifractal behaviors when we observe it at small spatio-temporal scales. This is a view point of the 
large-deviation scheme. 0j 

In the following we measure the distribution of the fractal dimension of the model on the square lattice at the critical 
temperature and show that the data is fitted by a non-Gaussian distribution derived from the Renyi or Tsallis entropy. 
Such a non-Gaussian behavior originates from the multifractality of the spin states and the q-index appearing in the 
distribution serves as a measure for the multifractality. On the other hand the non-Gaussian behavior is observed at 
small spatio-temporal scales as expected from the large-deviation scheme so that the g-index also serves as a measure 
for the degree of non-equilibrium. 

We measure the fractal dimension of the Ising model on the square lattice, 

LxL 



H = -JJ2°^> (!) 

<ij> 

where we adopt the ordinary notation, < ii > means the nearest neighbor bond and L is the number of the spins in 
x- or y-direction. The fractal dimension, || 

D = ln[M(l)/M(6)]/ln(6), (2) 

is determined by the Monte-Carlo data at the critical temperatute, T c = 2J/ln(l + ^/2)ks- Here b is the scaling 
parameter for the coarse-graining and LxL spins are transformed into (L/b) x (L/b) block spins. M(b) is the total 
spin at the level of b of the coarse-graining. The Monte-Carlo up-dating is done with the Metropolis algorithm. The 
simulation time, t, is represented by the unit of the usual Monte-Carlo step (MCS). In the following the measurement 
is done during r < t < r + T. The simulation is performed on our personal computer with Pentium-4 CPU. 

Some typical data for the unnormalized distribution of the fractal dimension are shown in Fig. 1 where D is measured 
at each Monte-Carlo steps and the total number of the measurement is T. The simulation data are well fitted by the 
distribution function, 

P(D) =P -[l + (q- 1)(D - D f/2a 2 ] l ^\ (3) 

which is derived by the maximum-entropy principle using the Renyi or Tsallis entropy (7[j^] so that the distribution 
exhibits a multifractality of the system. The Tsallis entropy is non-extensive, while the Renyi entropy is extensive. 
The functional form of the distribution is independent of the choice of the entropy, the Renyi or Tsallis, || so that our 
result is independent of the extensivity of the entropy. Here q represents the multifractality of the system investigated 
and the mean value, D ~ 1.875, is the one expected from the mono-fractal analysis. Q The limiting case with q — > 1 
corresponds to the Gaussian distribution which is parabolic in the D-\nP plane. The q-index for the solid curve is 
determined by the least-square fitting. In our scheme a is also determined by the data fitting. The peak value of the 
unnormalized distribution is Pq. 

In Fig. 2 the system-size dependence of the g-index is shown. 

In Fig. 3 the simulation-time dependence of the q-index is shown. 
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From the data in Figs. 2 and 3 we can read a tendency that q — > 1 at large spatio-temporal scales. The value, 
q = 1, is expected for the equilibrium state in the thermodynamic limit. Our data at small spatio-temporal scales 
reflect a deviation from this limit and exhibit non-equilibrium nature. A similar analysis has already been done for 
the experimantal data for turbulene. [Q 

In conclusion we numerically observed a multifractal behavior of the two-dimensional Ising model at small spatio- 
temporal scales. The g-index serves as a measure not only for multifractality but also non-equilibrium in this case. 
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FIG. 1. The distripution of the fractal dimension for a randomly-produced initial state with L = 125, 6 = 5 and r = 10°. 

FIG. 2. The system-size dependence of the q-index averaged over 10 randomly-produced initial states with b = 5, r = 10 5 
and T = 10 4 . 

FIG. 3. The simulation-time dependence of the q-index averaged over 10 randomly-produced initial states with L — 125, 
b = 5 and T = 10 5 . 



2 



I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 



2.0 - 



1.5 - 



m i 



1.0 



100 



200 



300 



400 



500 
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